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r-| ! Abstract : The nucleon structure functions are studied with a special emphasis on diquark 
correlations inside the nucleon. Assuming a phenomenological diquark-quark model of the 
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^ ■ nucleon at the low energy hadronic scale, we calculate the leading twist contributions to the 
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nucleon structure functions. We take the Nambu and Jona-Lasinio model for the calculation 
of the quark distributions in the diquarks, and use the convolution method for the diquark 
scattering process. The resulting quark distributions are evolved to the experimental mo- 
mentum scale using the perturbative QCD. We discuss the flavor dependence of the nucleon 
structure functions and the importance of the quark-quark correlations for the understanding 
of the nucleon structure. 
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1 Introduction 

Recent measurements of deep inelastic scattering show clear flavor dependences of the 
nucleon structure functions[l-3]. One example is the ratio of the neutron to proton structure 
functions F^(a;)/F2'(a;) which shows a large deviation from the naive quark-parton model 
prediction, which is 2/3. In addition, the recent NMC data shows that the difference between 
the dependence of -Ff (2^) ^'^'^ ^2^^) substantially larger than the perturbative QCD 
prediction H]. 

On the other hand, the diquark-quark model of the nucleon has attracted considerable 
interests. The idea of diquarks, i.e. correlated states of two quarks, is introduced phenomeno- 
logically to explain the scaling violation of the nucleon structure functions, which may be 
caused by the non-perturbative diquark correlations [Q. The diquark model also reproduces 
the proton-neutron differences of the structure functions]^]. Such a diquark structure is due 
to the quark-quark correlations which depend on its flavor structure, and may produce an 
asymmetry of the momentum distribution in the nucleon. In fact, such asymmetric momen- 
tum distributions are found in the calculations of the QCD sum rule|P, 0, where the most 
part of the proton momentum is carried by the u-quark with its spin directed along the 
proton spin and the remaining small part is carried by the u-d quarks with combined spin-0. 
Additionally, the elastic form factor of the nucleon requires an asymmetric momentum dis- 
tributions, which indicates the diquark clustering inside the nucleon This model is also 



supported by the study of several high energy scattering processes[0, 10 



In the low energy side of the hadron physics, the diquark-quark structure of baryons also 
plays various roles |[10|, 0, This picture was flrst introduced to explain the universality 



of the Regge trajectories|ll]. In particular, a recent work of Stech and his collaborators 
showed that the diquark correlations, which are the strongest in the scalar diquark channel, 
can account for the A/ = 1/2 problem of weak non-leptonic transitions, and reproduce the 
experimental data very well|T3|, |T^. Hence, it is of our great interest to study the flavor 
dependence of the nucleon structure in terms of the diquark-quark model as a successful 
low energy model of baryons. We note that such an approach naturally reproduces the 
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asymmetry of momentum distributions in the nucleon obtained by the QCD sum rule. 

Several theoretical calculations of the nucleon structure function were made in terms of 
the low energy quark models [15-20]. The structure function at the low energy scale = Q^, 
where the phenomenological quark model is supposed to work, is obtained by calculating the 
twist-2 matrix elements, since these matrix elements are dominant over higher twist terms 
in the Bjorken limit. The resulting structure functions are evolved to the experimental scale 
with the help of the perturbative QCD. Thus, one can compare them with experimental 
data. Here, we use the same method to evaluate quark distributions within the diquark 
model. 

Recently, the importance of the quark-diquark intermediate state in the calculation of 
the nucleon structure function was discussed within the MIT bag model ||19| and the quark- 
diquark modelpO|. It was pointed out that the mass difference between the scalar(O^) 
and the axial-vector(l'*') diquarks is crucial to account for the flavor dependence of quark 
distributions, since this mass difference produces the dominance of u-quark distribution in 
the proton at large Bjorken x. In such works, the diquarks are treated only as a spectator. 
However, if the two quark correlation is so strong as obtained in the QCD sum rule, we 
should take into account the diquark scattering process explicitly with the other quark being 
a spectator. 

In this paper, we shall calculate the nucleon structure function in the framework of the 
SU{6) diquark-quark model. Thus, both quark and diquark scattering processes contribute 
to the structure function, which are illustrated in Fig. 1(a) and (b). The process of Fig. 1(a), 
in which a quark is struck out by the virtual photon with the residual diquark being a 
spectator, is evaluated by using the standard method in the impulse approximation[2C]. 
For the diquark part, we use the convolution technique [^, ^ to calculate contributions of 
the diquark process shown in Fig. 1(b). We obtain the diquark terms as the convolution of 
quark distributions in diquarks with the diquark distributions in the nucleon. This diagram 
represents the leading twist contributions of diquarks, in which diquarks break up completely 
after absorbing the virtual photon^]. Hence, the resulting quark distributions scale; i.e. 
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there is no dependence without the QCD radiative correction. Here, we shall not calculate 
the twist-4 non-scaling contributions of diquarks, which are under consideration ||2^ . 



2. Nucleon Structure Function in the Diquark model 



The SU{Q) diquark-quark wave function of proton is written as [112 



\p t>= -^[3S{ud)u T + 2A{uu)+d j - V2A{uuyd t 
V 18 

- V2A{ud)+u i + A{udyu ]] (1) 

where S{ij) denotes the scalar (spin singlet; 0^) diquark with the flavor content i and j, and 
A{ij)'^ the axial-vector (spin triplet; 1"*") diquark with the helicity k. In this model, both 
the quark and the diquark scattering processes contribute to the hadronic tensor. Using (|T|), 
the proton and the neutron structure functions are given by. 



= l<l'ix) + ^g^(x) + ^QU^) + ^g^(x) (2-a) 
F,^x) = ^q'{x) + ^ix) + ^QUx) + IqUx) (2-b) 



where and are the quark distributions with the residual diquarks being the scalar and 
the axial-vector diquarks, respectively. Qfy and QY) are the quark distributions obtained 
by the scalar and the axial-vector diquark scattering processes, respectively. We can also 
express the valence distributions of the u and d quarks in the proton as. 



u{x) 



Iq'ix) + ^g^(x) + ^gg(x) + ^Qlix) 



d{x) = -q^ (x) + -QUx) + -Q 



D 



X] 



(3-a) 
(3-b) 
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Note that q^{x) term contributes only to the u-quark distribution. The flavor dependence of 
the structure functions arises from the difference of the distribution contents in this model. 
We also note that the definition ( p-a| , p-b|) , (|3-a|j3-b|) is different from that of the previous 
workp, in which the diquarks are assumed to be point-like. We consider the scaling 
contributions of diquark processes, where the diquarks are broken up completely by the 
virtual photon. 



3 Nambu and Jona-Lasinio model 

We use the Nambu and Jona-Lasinio (NJL) model to take into account the quark corre- 
lations, where the chiral invariance is the main ingredient [25]. Recently, the NJL model is 
studied in various subjects of hadron physics as a low energy effective theory of QCD [p6| . 
This model describes the SU{3) / meson properties very well with the parameters fixed by the 
pion and the kaon properties, in spite of the lack of confinement. The SU (3) NJL lagrangian 
for mesons is given by 26|, 



^NJL — + 

Co = ipiiYdf, - mo)^} 
Cm = Gfli^ta^y + i^tat^'^)^] 

-G^m^j,^f + i^Pta7,7'^PY] (4) 

Here ip denotes quark fields and mo their bare masses, ta are the flavor SU (3) / operators with 
the normalization tr{tatb) = Sabf^, and Gg^, Gy are the coupling constants. In this model, 
the quarks acquire the constituent masses dynamically due to the spontaneous breakdown 
of the chiral symmetry. The meson masses are obtained by solving the Bethe-Salpeter (BS) 
equation. 

On the other hand, the NJL lagrangian for the quark-quark interaction at the Fierz 
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transformed level is given by [26, 27], 

Cd = Gs[{i!taCAi'c){i'ctaCA^) + ii'taCAh^^c){i'ctaCAi'l^^)] 
- Gy[{ljjtsCA'j''i'c){i'ctsCA'lf,1p) 

+ i^taCAl''l^^c)i^ctaCAl^l^^)] (5) 

Here, ipc = Cip'^ are the charge conjugate quark fields and ca the antitriplet color SU{3) 
operators with the normalization tr{ciCj) = Sij/2, where A runs over only A = 2, 5 and 
7. The fiavor suffix a runs over a = 2, 5 and 7 for pseudoscalar, scalar and vector diquark 
terms, while it runs over s = 0, 1, 3, 4, 6 and 8 for the axial-vector diquark term due 
to the antisymmetrization of diquarks. Note that the presence of the charge conjugation 
operator produce an additional minus sign to the parity of diquarks. We get the masses 
and the wave functions of diquarks by solving the quark-quark BS equations as the case 
for mesons. Recently, the diquark structure is extensively studied in the framework of the 
NJL model||26|. We find that the structure of diquarks is quite similar with that of mesons 



with the corresponding quantum numbers. In particular, the scalar diquark [ipJl^'^) , which 
corresponds to the pion {tpij^ip) , shows a strong quark-quark correlation, whereas the axial- 
vector diquark {ipc'^^ip), which corresponds to the rho meson {ip'-j^ip), is a weakly bound 
quark-quark state. If the coupling constants satisfy the condition G^^ = Gg and Gy = Gy, 
the diquark mass is the same as the corresponding meson mass, which is so-called the Pauli- 
Giirsey symmetry between mesons and diquarks |]14|, |2^, Note that the strong diquark 



correlation in the scalar channel causes the large enhancement of the A J = 1/2 transition 
matrix elements of the hyperon non- leptonic weak decays |]T3|, |n . 



The nucleon properties such as the mass and the wave function within the NJL model 
should be obtained by solving the relativistic three body problems, and such an effort is being 
made with the relativistic Faddeev method p8[|. It seems, however, necessary to incorporate 
confinement for a successful description of baryons, which is absent in the NJL model. 
Therefore, we simply assume the diquark-quark model for the nucleon, and consider the case 
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of diquark-quark vertex being scalar pO|; 



(6) 



where 1 is the unit matrix in the Dirac space and (^(p^) is a regularization dependent function 
to be specified later. The diquark properties are obtained by solving the NJL model, which 
tells us the information of the quark-quark correlations. 



4. Calculations of Hadronic Tensor 

First, we evaluate contributions shown in Fig. 1(a), where a quark is struck out by the 
virtual photon with the residual diquark being a spectator. This part can be calculated 
following the work of Meyer and Mulders|2y]. We use the impulse approximation, and thus 



the hadronic tensor is represented by an incoherent sum of various processes. We define the 
constituent quark mass m and the diquark mass 1710 inside the nucleon, though the diquark 
and the quark are not the eigenstates of QCD. Their values are obtained within the NJL 
model. The calculation of the hadronic tensor in the Bjorken limit yields pO|, 



dp 

(27r)4(p2_^2)2^^' 

X -^—Tr\ 

dpi i^^M^ + 2mM-ml)+m^ + {l-x)pl] (7) 



X ^J^'^^ii^ + + m)(f + m)] 



2^2 I 



where 

Plmin = -r—rnjj - xM'\ 

I — X 

and M is the nucleon mass. Here, P is the proton momentum, and p and p2 are momenta 
of the struck quark and the spectator diquark, respectively. 

As for the diquark process shown in Fig. 1(b), we fold the quark distributions in diquarks 
with the diquark distributions in the nucleon. We use the NJL model for the calculation 
of the quark distributions in the diquarks. Note that these contributions to the structure 



functions scale. The approximation of the convolution method may be good for the scalar 
diquark process, since the scalar diquark is a small object due to its strong correlation. In 
fact, the NJL model calculation gives a small radius < 0.4fm. Phenomenological estimates 
of the scalar diquark radius also provides small value 0. But the size of the axial- vector 
diquark is comparable with that of the nucleon, thus the convolution is not suitable for the 
calculation of the axial-vector diquark process. Assuming the quark-diquark vertex to be 



scalar as mentioned above, we obtain the diquark term as the convolution integral[22 



Q'd{x) = / dx/y Fo/Niy) Df)i{x/y) (8) 

Jx 

Fo/Niy) is a probability to find a diquark in the nucleon with the light-cone momentum 
fraction y. This part is evaluated in terms of the scalar vertex (6); 



^D/Niy) = yjT^.^^k! r^2_^2,2 



X 2'K6{kl - m^)6{y - '—)e{k+)e{M - k+) 



Jtmin IGtC^ [t + mj^y 



where 

tmin = —^rri^ - yJVP. 

Here, the diquark carries the momentum q and the residual quark q2. D£,i{x) is the quark 
distribution of the leading twist diquark structure function. For the calculation of Dciix), we 
use the same procedure as done in the meson case within the NJL model [29]. For example, 
the quark distributions in the scalar diquark is given by, 



2 

Dsix) oc gl^^ dfil[ +x /".J fif^l 



X 6'(m|,x(l - x) - xm^ + (1 - . (10) 
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gsqq is the diquark-quark-quark coupling constant obtained by the NJL model calculation. 
fifi'E) is the regularization function. Note that the form of (|1^) is the same as the quark 



distribution of the pion obtained in ref . [£9| , since the diquark structure is the same as that 
of the corresponding meson in the NJL model due to the relation p6|, p7| ; 

where C is the charge conjugation operator and Sp the quark propagator. The structure 



function of the axial- vector diquark is also calculated in the similar manner pO| . Thus, we 
get the convolution of the axial-vector diquark term. 

The integrals (0), (P), and (|l^) go to infinity due to the non-renormalizability of the 
NJL model. To avoid the divergence, we introduce the regularization function /(p^) in the 
Euclidean space. There exists an difficulty in the introduction of the standard NJL cutoff 



for DIS process |pl|, |3^ . Here, we use the Fermi distribution type cut-off function, which is 



consistent to the usual sharp cut-off method except for large momentum. 

Here, we use Ar^lGeV and a ~ O.lGeV^ to reproduce the meson properties. We examine also 
other plausible form like the exponential cut-off for completeness. We obtain similar quark 
distributions except for the shape around x ~ 1, which depends on the high momentum 
behavior of the cut-off scheme. 



5. Results and Discussions 

We present here the calculated result on the structure function. Concerning the NJL 
model parameters, we take the current quark mass rrio = 5.5MeV as a semi-empirical value, 
and GfA"^ = 2G^'A'^ = 22.2, A = 860MeV, a = 0.15GeV^ which are fixed by the pion 
properties. For the diquark part, we treat the diquark masses as free parameters, and study 
their dependence of the structure function. By the analysis of the N-A mass splitting in the 



one-gluon exchange picture, the scalar diquark mass ms was assumed to be 575MeV"[33 
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and the axial- vector diquark mass to be 200MeV higher than the scalar di quark mass|T9|, 
20| . Much larger scalar- vector diquark mass difference was obtained in the instanton liquid 



modell^l] or other models|T^, Here, we assume = 775 Me V", which is consistent to 
the above arguments and comparable with the vector meson mass like p. Then, we investigate 
the scalar diquark mass dependence of the results. The quark correlation affects the quark 
distributions through the diquark masses. 

We first show the quark distributions of various processes at the low energy scale in Fig.2. 
These distributions represent their strength of the correlation. The quark distribution q''^{x) 
(solid curve), where the scalar diquark is a spectator, has a sharp distribution peaked at 
X ~ 0.5, because of the weakness of the interaction between the scalar diquark and the third 
quark due to the strong binding of the diquark in the scalar channel. q^{x) (dashed curve) 
shows slightly a broader distribution. q^{x) and q^{x) are almost unchanged from the results 
of ref. [^. In addition, we have calculated newly the contributions from diquark process 
(Fig. 1(b)). The momentum of two quarks in the scalar diquark is, compared to q''^{x), peaked 
at small x; a; ~ 0.25, as shown by dotted line. Using these distributions, we can obtain the 
momentum ratio carried by each quark. For the scalar channel, we find at 7715 = 575 Me^ 

< xq^ >:< xQf) >~ 2.2 : 1 , 

where < xq >= Jq dxxq{x). This means that the u-quark with the scalar diquark being a 
spectator carries the largest part of the proton momentum, and the remaining small part is 
carried by the u-d scalar diquark due to the strong correlation in the scalar channel. This 
result is consistent with the QCD sum rule[§, 0. If the diquark correlation becomes weak, we 
find that this ratio is close to 1. For example, < xq^ >:< xQf) >~ 1.3 : 1 at ms = 775MeV. 
In the case of axial-vector diquark process, we obtain 

< xq^ >:< xQY) >~ 1.3 : 1. 

Thus, both two quarks in the axial- vector diquark {Q^{x): dash-dotted curve) and the third 
quark {qv{x)) have almost the same amount of the proton momentum. This fact reflects 
the weak correlation of quarks in the axial- vector diquark channel. Comparing Qf){x) with 
Q^{x), the difference of both distributions is not so evident. Originally, the shapes of 
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distributions in the scalar and axial- vector diquarks are quite different, as tt and p in the 



meson case 



l30[| . However, such a difference is washed out after the convolution with the 



entire diquark motion in the proton. 

We take the low energy model scale Ql = 0.2GeV^, which is used in ref. |^. At this 
scale, the running coupling constant is still small; as{Ql)/n ~ 0.4. Thus, the inclusion of the 
second order QCD corrections gives a small change for the evolution from our result within 
10%[Q. We use the first order Altarelli-Parisi equation [p6|] with the Aqcd = 250MeV, in 
order to compare our result with experiment. For the sea quark and the gluon parts, we 
simply assume that these distributions vanish at Q'^ = Ql, and are generated by the QCD 
evolution process. 

The proton and the neutron structure functions -Ff (^) -F^(x) thus obtained at = 
15GeV'^ are shown in Fig.3 with experimental data|jl|, where we use rriA — ms = 200MeV. 
Our result shows a reasonable agreement with data for x > 0.4, and is somewhat different 
for small x value. One of the possible solutions for this disagreement is to take into account 
the sea quark degrees of freedom explicitly at the low energy scale. The inclusion of the 
q — 4q or q — 3qq intermediate state at the model scale may change the shape of F2{x) 
at small x as introduced in ref. ll^. Comparing our calculated results with the ones of 



the MIT bag model[^ or the diquark spectator model||2^, our structure functions are 
distributed over much wider range of the Bjorken x. This fact is the consequence of the strong 
diquark correlation. We also show in Fig.3 the proton structure function without the diquark 
scattering part (Fig.l(a) only), in which we use Ql = 0.2GeV^. This result is essentially 
the same as that of the previous model calculation without the diquark process[16-20]. This 
structure function has a sharp peak around x ~ 0.3, and disagrees with experiment. In order 
to reproduce the absolute values and the peak position of the experimental data, one must 
take a smaller value for the model scale Q^. In fact, a much smaller value for the low energy 



scale is used in the previous calculations |T7|, [18|, ^ so as to obtain a better agreement with 
experiment. 

We would like to discuss the dependence of the structure function on the low energy scale 
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Ql- The nucleoli structure functions are shown in Fig. 4, with the choice of the low energy 
scale = 0.125GeV^. The running coupling constant at this scale is about two times 
larger than that for Ql = 0.2GeV^. As compared with the result in Fig.3, the calculated 
structure functions are in a good agreement with experiment, especially at low x. However, 
the use of the perturbative QCD becomes less reliable at such a low Q^. 

We discuss how the flavor structure of the structure function depends on the diquark 
correlation. The ratio F2{x) / F2{x) is shown in Fig. 5. Our result is in good agreement with 
experiment ^ |^, in the case niA — rris = 200MeV. This is due to the dominance of u- 
quark distributions, namely q^{x), at large x, which is caused by the asymmetric momentum 
distributions of quarks and diquarks. In the middle x range, x ~ 0.4, the resulting ratio is 
close to 2/3, since this part of distributions is dominated by the diquark parts, i.e. Qf){x) and 
QY){x). We find that our calculated result is somewhat smaller than experiment for small x. 
The inclusion of sea quark at the low energy scale enhances this ratio in the small x region, 
and may resolve this discrepancy. We also show in Fig. 5 the case mA—Tns = 0, 100, 275MeV. 
If we take the same values for the scalar and the axial-vector diquark masses, the ratio is 
close to 2/3. We also show the difference -Ff (^) ~ F2{x) in Fig. 6. The peak position of our 
result is consistent with experiment, but the absolute value is rather large. This behavior 
was found in the previous calculations in terms of other models. 

6. Conclusion 

In conclusion, we have studied the flavor structure of quark distributions based on the 
phenomenological diquark model. This model reproduces the asymmetric momentum dis- 
tributions in the nucleon obtained by the QCD sum rule. The neutron to proton structure 
function ratio F2{x)/F2{x) shows a good agreement with the available DIS data. The agree- 
ment is lost as the diquark correlation becomes weak. Our result indicates that the quark 
correlation, especially in the scalar channel, is important to understand the DIS data, which 
is also crucial for the low energy hadron properties such as the AI = 1/2 non-leptonic weak 
transition. Although we consider the twist-2 contributions of the structure function in this 
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paper, twist-4 pieces of them give us much information about the quark-quark correlation |p7||. 
The diquark-hke correlation may be important to study such higher twist effects. Work along 
this line is in progress. 
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Figure Caption 

Fig. 1 : The forward scattering amplitude ("liandbag diagram") of the nucleon. The thin 
solid line represent the quark, and the shaded line the diquark. The nucleon and the virtual 
photon are depicted by the thick solid and wavy lines. For details and notation, see text. 

Fig. 2 : The quark distributions of the nucleon at the low energy scale = Ql as 
a function of the Bjorken x. The solid and dashed curves denote the quark distributions 
obtained by the calculation of Fig. 1(a) with the residual diquarks being scalar or axial- vector, 
respectively. The dotted and dash-dotted curves represent the distributions of the scalar and 
axial- vector diquark scattering processes (Fig. 1(b)). 

Fig. 3 : The proton and neutron structure functions at = IbGeV"^. The solid curve 
represents F2{x), and the dashed curve -F^(x). The proton structure function calculated 
without diquark scattering processes is depicted by the dotted curve. Here, we use the 
low energy model scale = 0.2GeV^. The experimental data are taken from the EMC 
experiment JI]] . 

Fig. 4 : The proton and neutron structure functions at = 15GeV^. The notations are 
the same as those of Fig. 3. The low energy scale Ql = 0.1256*6^^ is used. 

Fig. 5 : The ratio of the nucleon structure functions F2/F2 at = 15GeV'^ with 
mA-ms = 0, 100, 200, 275 MeV^ {Ql = 0.2^6^^). The experimental data are taken from 
refs. [1-3]. 

Fig. 6 : The difference of the nucleon structure functions at = IbGeV"^ with mA—^s = 
0, 100, 200, 275MeV {Ql = 0.2GeV^). The experimental values with error bars are taken 
from the EMC and BCDMS experiments [|l|, ||. 
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